Activation cross-sections of longer lived products of deuteron induced nuclear 

reactions on ytterbium up to 40 MeV 



F. Tarkanyi'', F. Ditroi" *, S. Takacs'', A. Hermanne'', H. Yamazaki'', M. Baba% A. Mohammadi'', A.V. Ignatyuk'' 

"Institute of Nuclear Research of the Hungarian Academy of Sciences (ATOMKI), Debrecen, Hungary 
^ Cyclotron Laboratory, Vrije Universiteit Brussel (VUB), Brussels, Belgium 
Cyclotron Radioisotope Center (CYRIC), Tohoku University, Sendai, Japan 
'^Institute of Physics and Power Engineering (IPPE), Obninsk, Russia 



Abstract 

In the frame of a systematic study of the activation cross-sections of the deuteron induced nuclear reactions, excitation 
functions of the "«'Yb(d,xn)i^'''i^3,i72mg,i7img,i70,i69L^^ ™'Yb(d,x)i^5'i^''Yb and ««'Yb(d,x)i''3'i^2,i68,i67,i65rpjj^ reactions 

are studied up to 40 MeV, a few of them for the first time. Cross-sections were measured with the activation method 
using a stacked foil irradiation technique and high resolution y-ray spectrometry. The experimental data are analyzed 
and compared to the results of the theoretical model codes ALICE-IPPE, EMPIRE-II and TALYS. From the measured 
cross-section data integral production yields were calculated. Applications of the new cross-sections are discussed. 
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1. Introduction 

Activation cross-sections data of deuteron induced 
nuclear reactions are important both for application in 
several fields and for benchmarking description of re- 
actions with different model codes. In connection with 
several projects it was recognized that a high priority 
needs to be given to establish a reliable database for 
deuteron activation data. Literature search shows that 
the status of the available cross-section data for deuteron 
induced reactions (especially above 15-20 MeV) is very 
poor No systematical study has been done earlier and 
the cross-section data (except a few well measured mon- 
itor and medically important reactions) show large dis- 
crepancies. To meet the requirement of improving the 
reliability of available data, we started to establish an 
experimental activation database some years ago, by 
performing new experiments and a systematical survey 
of published deuteron induced activation cross-sections 
up to 50 MeV. The targets were irradiated with external 
beams of the Debrecen, Brussels, Louvain-la-Neuve, 
and Sendai cyclotrons. The mostly new, measured exci- 
tation functions are compared with the results of differ- 
ent nuclear reaction model codes. As a continuation of 



our investigation of deuteron induced activation cross- 
sections on rare earth targets we performed an irradi- 
ation of Yb targets with a 40 MeV incident deuteron 
beam at the Sendai cyclotron. Due to the experimental 
circumstances (rather long cooling times before start of 
activity measurement) this experiment results in excita- 
tion functions for longer lived Yb, Lu and Tm radionu- 
clides. We have already earlier experimentally studied 
most of these activation products in proton and alpha- 
particle induced nuclear reactions on ytterbium lUQ • 

2. Earlier investigations 

Only a few earlier experimental works have been 
found dealing with deuteron induced activation products 
on Yb, mostly at lower energies. 

• Nichols et al. investigated the excitation functions 
of Yb(d,xn) Lu, i^-Lu, '^^Lu between 1 3.4 to 29.0 
MeV a. 

• A. Hermanne et al investigated the cross sections 
for deuteron-induced reactions on Yb and Mea- 
sured cross-sections between 3 and 20 MeV for 
Yb(d,xn)i^"Lu/'^iLu/i^2Lu/i"Lu/i^'*Lu/'^^Lu, and 
Yb(d,xnp)i''''Yb/i^5Yb 0. 
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• S. Manenti et al. measured the 

activation cross-sections of Yb(d, 
xn)i6''Lu/i™Lu/i^iLu/i^2Lu/i73Lu/i74Lu/i76Lu/i77Lu^ 

Yb(d,xnp)i®Yb/i''5Yb/i^^Yb reactions up to 18.18 
MeV |6|. 

• P.P. Dmitriev et al., in the frame of a systematic 
study: "Yields of Radioactive Nuclides Formed 
by Bombardment of a Thick Target with 22-MeV 
Deuterons" measured the thick target yields on Yb 
for production of '^^''^"^Lu |7|. Considering the 
theoretical results on production cross-sections for 
deuteron induced reactions on different target ma- 
terial we have found also only few works: 

• Results calculated by the ALICE-IPPE code M 
can be found in |5|. 

• The most recent results of systematic calculations 
made by the TALYS code (91 can be found in the 
TENDL 2012 on line library 001 . 

• Comparison of evaluated and experimental data 
from ENDF/B-Vll.l, TENDL-2011 and EXFOR 
can be found in the JANIS Book of deuteron- 
induced cross-sections prepared by N. Soppera et 

aiim. 

3. Experimental 

Elemental cross-sections on '""Yb targets were mea- 
sured using the activation method and the standard 
stacked foil irradiation technique combined with high- 
resolution y-ray spectrometry. The irradiation, the ac- 
tivity measurement and the data evaluation were similar 
as described in more detail in our recent works lfT2ll . 
The excitation functions were measured up to 40 MeV 
by bombarding ""'Yh targets with low intensity deuteron 
beams at the AVF-930 cyclotron of the Tohoku Univer- 
sity (Sendai). Reactions induced on Al foils present in 
the stack were used to monitor the parameters of the 
bombarding deuteron beams. The ytterbium (thickness 
23 fim) and aluminum foils (thickness 100 fim) were of 
high purity ( 99.5%) and purchased from Goodfellow 
(UK). The irradiated stack had a complex structure, con- 
sisting of a 9 times repeated sequence of Rh(12.3 fJ.m), 
Ho(25 fim), Au(10.7 jum), Yb(23 fim ), MoRe alloy(50 
yum), CuMnNi alloy(25 fim) foils and the Al(100 /urn) 
foils serving as monitor. In principle (n,y) and (n,2n) 
reactions induced by secondary neutrons can also con- 
tribute to the production of the isotopes of Yb. We have 
checked the possible effect of neutron induced reactions 
by inserting additional target foils at the end of the stack 



holder, where complete stopping of the deuterons is ex- 
pected. In the used irradiation setup the effect of neu- 
trons was found to be negligibly small. The results for 
the Mn, Rh, Au and MoRe targets were already pub- 
lished fTJ-TSl and for Ho and Ni target material will 
be discussed in separate reports. The irradiation of the 
targets was performed in He-gas atmosphere in a water 
cooled target holder For determination of the beam in- 
tensity and energy, the complete excitation functions of 
the ^^Al(d,x)^^'^'*Na monitor reactions were measured, 
simultaneously with the excitation functions of reac- 
tions induced on Yb targets (see Fig. 1). The target 
stack was irradiated with a coUimated deuteron beam of 
40 MeV incident energy, for 30 min at about 24 nA. 




Energy (MeV) 



Figure 1: Application of monitoring using the Al(d,x) • Na reac- 
tion 

The radioactivity of each sample and monitor foil was 
assessed nondestructively by HPGe y-spectrometry in 
two series of measurements. The first counting series 
was started about 20 hours after the end of the bom- 
bardment (FOB) (each foil was measured for 10 min), 
the second series started about 20 days after EOB (each 
foil was measured for 5-8 h). The large number of ir- 
radiated foils and the high induced activity prevented 
the optimization of the spectrum measurements of the 
samples. The information on activity of short-lived ra- 
dionuclides was lost due to initial one day waiting time 
after EOB. Due to the large number of foils and the lim- 
ited number of detectors (2), the measuring time of the 
first assessment was restricted to 10 min. An additional 
inconvenience of these circumstances is that detection 
of the weak y-lines of medium half-life radioisotopes 
was also impossible. The decay data were taken from 
NUDAT 2.6 [16 1, the standard cross-section data of 
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the used monitor reactions ""'Al(p,x) • Na from IAEA 
recommended database |T7|. For the reaction Q-values 
we used the NNDC Q-value calculator [18J. The energy 
degradation along the stack was determined via calcu- 
lation 1 19 1 and corrected on the basis of the simulta- 
neously measured monitor reactions by the method de- 
scribed in ll20l . The excitation functions were hence 
determined in an accurate way relative to the simulta- 
neously re-measured monitor reactions. The decay data 
and the contributing reactions with their Q-values are 
presented in Table 1. The naturally occurring ytter- 
bium is composed of seven stable isotopes, therefore, 
so called elemental cross-sections were determined ex- 
cept for '^^Lu, where only (d,p) and (d,n) reactions on 
'^^Yb can contribute (see Table 1). The uncertainties on 
the cross-section values were estimated in the standard 
way: the independent relative errors of the linearly con- 
tributing processes were summed quadratically and the 
square root of the sum was taken ETl . The uncertainty 
of the energy of the bombarding beam in the middle 
of each irradiated foil was estimated combining the un- 
certainties of the energy degradation depending on the 
thickness of the target foils in the stack and the energy 
straggling. 

4. Comparison with the results of model codes 

The experimental data are compared with the cross- 
section data reported in the last two TALYS based 
TENDL 2011 and TENDL 2012 Activation Data Li- 
braries IflOl to show the agreement with the experi- 
mental data and contribute to the development of the 
TENDL data library. The cross-sections of the inves- 
tigated reactions were calculated by us,using the mod- 
ified pre-compound model codes ALICE-IPPE |8| and 
EMPIRE-II |22|. hi the modified ALICE IPPE-D and 
EMPIRE -D code versions the direct (d,p) channel is in- 
creased strongly, which gives better agreement for the 
(d,p) channel, but naturally it reflects also in the results 
for all other reactions Il23ll24l . 

5. Results 

The measured experimental cross-section data are 
shown in Figs 2-15 together with the earlier literature 
results and results of the theoretical calculations. The 
numerical values are collected in Tables 2-3. The fol- 
lowing notations were used: 

• the "mg" denotes the activation cross-section of the 
ground state after the "complete" decay of the sig- 
nificantly shorter half-life isomeric state decaying 
partly or completely to that ground state; 



• the "cum" denotes the activation cross-section of 
the final product after the "complete" decay of the 
simultaneously produced significantly shorter half- 
life parent nuclei decaying partly or completely to 
the investigated final product. 

5.1. Cross sections oflutetium radioisotopes 

Except for production of '^^Lu only direct (d,xn) 
reactions contribute to the formation of the measured 
lutetium radioisotopes. 

5.1.1. ^^^Yb(d,xn)^'^^Lu 

In this study, only the cumulative production of 
'^^'•'Lu (6.71 d half-life) following total decay of parent 
'^^Yb (Tl/2 = 1.9 h) was assessed. It practically con- 
tains no contribution from the internal decay of '^^'"Lu 
(160.4 d, IT 21.7%) (long-lived, low formation cross- 
section). The medically relevant '^^*Lu can be pro- 
duced in a direct reaction via '^^Yb(d,n)'^^*Lu or in- 
directly by decay of parent '^^Yb obtained through the 
''"'Yb(d,p)'''^"'^^Yb route. According to the earlier low- 
energy measurements ||5] [6| and to the results of the 
model calculations, the (d,p) reaction has a significant 
predominance over the direct (d,n) process. The cross- 
section data above 20 MeV obtained in this study, have 
large uncertainties due to the separation process of the 
208 keV y-line of '^^*Lu overlapping with the y-line 
at nearly the same energy of "'^Tm (with a half-life of 
9.25 d, similar to '^^"^Lu). Our new data are higher than 
the literature experimental data in the overlapping en- 
ergy region (Fig. 2). The TENDL 2011 values present 
total cross-sections of '^^Lu but in the TENDL 2012 al- 
ready the cumulative cross-section of the '^^-^Lu can be 
deduced. Taking into account the small cross-sections 
of the direct production of the ^"'"Lu the TENDL 201 1 
results reproduce well the experimental data, both the 
shape and the magnitude. In case of the TENDL 2012 
the underestimation of the experimental data is signifi- 
cant. 

5.7.2. ""'Yb(d,xn)^^^Lu 

The production of '^■'Lu (1.37 a) arises from reac- 
tions on four stable Yb isotopes. The contribution of 
the '^^Yb(d,n) reaction is small and cannot be distin- 
guished on Fig. 3. The i"Yb(d,2n) and ■^4Yb(d,3n) 
produce the first maximum and the '^*Yb(d,5n) reac- 
tion the second one (Fig. 3). The agreement with the 
earlier experimental data of Hermanne et al. [5 1 and 
Manenti et al. |6 1 is acceptable good, but there is signif- 
icant disagreement with the results of Nichols et al. ID. 
It is difficult to judge the agreement with the TENDL 
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Figure 2: Excitation function of tlie Yb(d,x) Lu process 

predictions: at the low energy maximum it seems to be 
good, but at higher energies the experimental data are 
higher ALICE-D and EMPIRE-D runs together with 
the TALYS results up to 15 MeV, at the maximum both 
give higher values, and at higher energies give similar 
results again. 



arises from contribution of the (d,4n) reaction on high 
abundance '^"^Yb. Our new experimental data show 
good agreement with the earlier results and give new 
data above 30 MeV. The both TENDL curves run al- 
most together giving good description of the experimen- 
tal results in the studied energy domain, while ALICE- 
D shows a larger maximum and EMPIRE-D overesti- 
mates above 15 MeV. 
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Figure 4: Excitation function of the '""YbCd.xn) Lu process 
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Figure 3: Excitation function of the '""YbCd.xn) Lu process 

5.7.5. ""'Yb(d,xn)™Lu 

The acceptable agreement of the experimental and 
theoretical excitation functions for '^^Lu production is 
shown in Fig. 4. Reactions on five stable Yb isotopes 
are contributing (see Table 1). The cross-sections con- 
tain the complete contribution of the decay of short- 
lived isomeric state (3.7 min). The large maximum 



5.1.4. "'"Yb(d,xn)^^^Lu 

Below 30 MeV a good agreement of the available ex- 
perimental data and the theoretical results was obtained 
(see Fig. 5) for cumulative production of '^'Lu (8.24 
d) (including complete decay through isomeric transi- 
tion of the short-lived (79 s) isomeric state), but our val- 
ues are systematically higher at higher energies. The 
three maxima can be connected to the reactions on re- 
spectively the low mass ytterbium isotopes, to the '^"^Yb 
and to the '^^Yb target isotopes. Our new data agree 
well with those of Hermanne |i5J and Manenti |6| up 
to 20 MeV and significantly larger than Nichols' re- 
sults above 20 MeV. The local maxima era also repro- 
duced well. Now the best approximation is given by 
the ALICE-D code. The TENDL results are also ac- 
ceptable under 20 MeV. EMPIRE-D overestimates the 
experimental data in the whole energy region. 

5.7.5. "'"Yb(d,xn)™Lu 

The excitation function for production of '^*'Lu 
(2.012 d) is shown in Fig. 6. '™Lu has no isomeric 
state, so the presented results are direct cross-sections 
of the ground-state, resulted in from (d,xn) reactions 
on different Yb stable isotopes (M=170-174). Our new 
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Figure 5: Excitation function of the '""YbCd.xn) Lu process 



experimental data (Fig. 7), while both EMPIRE-D and 
ALICE-D overestimate significantly above 20 MeV. 




Deuteron energy (MeV) 



data agree acceptable well with the result of earlier ex- 
periments below 20 MeV. The excitation functions pre- 
dicted in the TENDL show also a good agreement with 
the experiment. It is difficult to distinguish between the 
results of the different model calculations, while they 
mostly run together. Only above 30 MeV is obvious 
that both TENDL data underestimate the experimental 
results. 



Figure 7: Excitation function of the "'"Yb(d,xn)'*''Lu process 

5.2. Cross-sections of ytterbium radioisotopes 

The radioisotopes of Yb can be formed through two 
routes: directly via (d,pxn) reactions and by decay of 
simultaneously produced parent radioisotopes of Lu and 
Tm with the same mass. 
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5.2.7. ^'^''Yb(d,pf'^Yb 

Due to the long waiting time (around 20 h) after EOB 
before the first measurement we could not identify the 
y-lines of the '^^Yb (1.911 h) in the spectra. To illus- 
trate the poor predictivity of the model codes for the 
(d,p) reactions, in Fig. 8 we have compared the ex- 
perimental data of lis] |6l with the TALYS results in the 
TENDL libraries. The disagreement in the maximum 
value is around a factor of two. The results of model 
calculations still significantly underestimate the exper- 
imental data, although new improvements for descrip- 
tion of (d,p) reactions were implemented in the TALYS 
code. Phenomenological systematics of the (d,p) cross- 
sections and the problems of the theoretical descriptions 
are discussed in more detail in I 



Figure 6: Excitation function of the "°'Yb(d,xn) Lu process 



5.7.6. ""'Yb(d,xn)^^'^Lu 

No earlier experimental data were found for produc- 
tion of '«^Lu (32.018 d), probably due to the high effec- 
tive reaction threshold not covered in the energy range 
in ||5]|6l. The TENDL results reproduce excellently the 



5.2.2. ""'Yb(d,x)^^^Yb 

Fig. 9 shows the cross-sections for cumulative pro- 
duction of '^^Yb (4.185 d) via direct (d,pxn) reactions 
and from decay of '^^Tm (15.2 min) produced by the 
'^''Yb(d,2pn) process. According to the theory and to 
the systematics the (d,2pn) contribution is small. There- 
fore the cumulative cross-sections mostly reflect the di- 
rect (d,pxn) cross-sections. Our present data above 20 
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Figure 8: Excitation function of the Yb(d,p) Yb process 

MeV are new, in the overlapping region show agree- 
ment with the data of Hermanne |5 |. The earlier results 
of Manenti (S) are higher. Similar to the previous case 
the theoretical results strongly underestimate these pro- 
cesses, the form of the excitation function is reproduced 
by TALYS only. 



TALYS, ALICE-D and EMPIRE-D results is acceptable 
good, although it is difficult to distinguish between the 
theoretical codes, while they run together up to 25 MeV. 
Above this energy our new experimental data are better 
described with ALICE-D and EMPIRE-D. The earlier 
experimental data of Manenti [6J and Hermanne |i5J are 
in good agreement with our new results in the overlap- 
ping energy domain. 
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Figure 10: Excitation function of the ™'Yb(d,x)'*'''Yb process 
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Figure 9: Excitation function of the "'"Yb(d,x) Yb process 

5.2.3. "'"Yb(d,x)^^'^Yb 

The isotope i®Yb((Ti/2 = 32.018 d) is obtained 
through the direct production via (d,pxn) reactions and 
is also populated by decay of the simultaneously pro- 
duced parent '^^Lu (34.06 h). The measured cumula- 
tive activation cross-sections are shown in Fig. 10. The 
agreement with the earlier experimental data and with 



5.3. Cross-sections of residual radioisotopes of thulium 

The radioisotopes of thulium are produced directly 
via (d,2pxn) reactions, and/or through decay of simulta- 
neously formed parent Lu and Yb radionuclides. For di- 
rect production at low energies the a-particle emission 
dominates the possible (p,2pxn) (x> 2) reaction chan- 
nels. At higher energies the individual particle emission 
is more significant. The small contributions from the 
yS -decay of erbium radioisotopes can be neglected due 
to the low probabiUty of the (d,3pxn) processes. 

5.5.7. ""'Yb(d,x)^^^Tm 

According to the Fig. 1 1 the agreement of the exper- 
imental and theoretical data for the directly produced 
'^■'Tm (8.24 h) is not too bad, especially under 30 MeV, 
taking into account the complex particle emissions for 
production of this radionuclide. 

5.3.2. ""'Yb(d,x)^^^Tm 

The measured excitation function for production of 
'^^Tm (93.1 d) is shown in Fig. 12. The theory 
significantly underestimates the experimental values. 
The both TENDL calculations run together, there is 
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Figure 11: Excitation function of tiie "°'Yb(d,x) '^Tm process 

no improvement between the 2011 and 2012 versions. 
EMPIRE-D and ALICE-D give similar, but still under- 
estimating results. 



20 




I ■ ' -"^ ^ ' ' ' 1 

10 20 30 40 so 60 

Deuteron energy (MeV) 



Figure 12: Excitation function of the "°'Yb(d,x) °°Tm process 

5.3.3. "'"Yb(d,x)^^'^Tm 

The radionuclide '^^Tm (9.25 d) is produced directly 
and from the decay of short-lived parent ^^'Yb (17.5 
min). According to Table 1 and to Fig. 13, the thresh- 
old for production of the '^^Tm is high in the investi- 
gated energy range, but the cross-sections are small. We 
present cross-sections for cumulative production, after 
complete decay of '^^Lu. The resulted cross-sections 
contain large imcertainties, especially near the effective 
threshold, due to the need of correction for interfer- 



ence in the y-ray signal of ^°'Tm. The two measurable 
y-lines of '^^Tm have energies of 207.801(42 %) and 
531.54 keV (1.61%). As it was mentioned earlier the 
207.8 keV line overlaps with the 208.3662 keV y-line 
of '^^^Lu, which is experimentally practically insepara- 
ble because of the energy resolution of the used detec- 
tor. Therefore we have tried to use the 532 keV y-line of 
'^^Tm with limited success (low count rates, large scat- 
tering). Therefore we have separated the contribution of 
^'^*Lu in the y-spectra measured 20 day after BOB. It 
was possible due to the lower y-abundance of the 208 
keV y-line of '^^Lu, due to the shorter half-life and to 
the relatively low cross-sections for cumulative produc- 
tion of '^^*Lu on natural Yb. The experimental exci- 
tation function for "'"Yb(d,xn)'^^Lu(cum) (discussed in 
section 5.1.1) is also shown in Fig. 13. It is difficult 
to distinguish, which model code gives the better re- 
sults, but at least at lower energies under 35 MeV the 
approximations of EMPIRE-D and ALICE-D are bet- 
ter, but above this energy the experimental values reach 
the TENDL 2012 curve. 




Deutervn energy (MeV) 



Figure 13: Excitation function of the ""'Yb(d,x)'°'Tm process 

5.3.4. "<"Yb(d,x)^^^Tm 

The '^^Tm (30.06 h) is produced both directly 
through "'"Yb(d,pxn) reactions and indirectly from de- 
cay of short lived parent ^^^Yb (9.9 min). The effec- 
tive threshold is close to our maximum energy. The few 
measured cross-section points are shown in Fig. 14. 

6. Integral yields 

On the basis of our new experimental cross- 
section data and the results of the model calcu- 
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Figure 14: Excitation function of the '""Yb(d,x)'^^Tm process 



Figure 16: Integral yields for production of Yb and Tm radioisotopes 



lations we determined integral yields for produc- 
tion of the investigated reaction products. The re- 
sults are so called physical yields ll27l correspond- 
ing to instantaneous short irradiation. The calcu- 
lated yields for ««'Yb(d,xn)i"-i^^^'i^2mg''^i'"«'i™'i^''Lu, 
"'"Yb(d,x)i"'i®Yb and ««'Yb(d,x)i^3,i72,i68,i67-fj^ ^^.g 

shown in Fig. 15-16 in comparison with the few ex- 
perimental integral thick target yield data found in the 
literature. Dmitriev et al. reported radioactive thick tar- 
get nuclide yields at 22 MeV incident deuteron energy 
Q for the ™'Yb(d,x)i"Lu and '""Yb(d,x)i^4,?L^ j.g^^. 
tions. Our calculated yield curve for '^^Lu runs above 
the data point of Dmitriev at 22 MeV. 




Figure 15: Integral yields for production of Lu radioisotopes 



7. Summary and conclusion 

Twelve excitation functions of deuteron induced nu- 
clear reactions on Yb were measured up to 40 MeV. 
All data above 20 MeV, except for production of 
I7i,i72,i73j_^^ were measured for the first time. The 
experimental data were compared with the results of 
EMPIRE-D and ALICE-D calculations and data in the 
TENDL 2011 and 2012 libraries based on TALYS re- 
sults. The theoretical calculations describe the shape 
and the absolute values of the experimental data only 
with moderate success, and the (d,p) reaction is espe- 
cially significantly underestimated by the theory. Fur- 
ther improvement of the used model codes is required. 
The comparison with the low energy earlier experimen- 
tal data shows excellent agreement. The obtained exper- 
imental data provide a basis for improved model calcu- 
lations and for different applications in the field of nu- 
clear medicine and for estimating production yields of 
long-lived products for dose estimations and waste han- 
dling. Presently the ytterbium has very few industrial 
uses, but the applications are growing (in laser technol- 
ogy, fiber optic technologies, etc.). Among the radioac- 
tive nuclear reaction products '^^Lu II28H31 1, '^^Lu |33- 
l34l . '^''Yb f35|, i"Yb |36| and ""^Tm f371 have poten- 
tial relevance in nuclear medicine for cancer treatment 
and as tracer in nuclear biology (for a review see 138]). 
In the domain of industrial applications, the y-emitter 
'^^ Yb (39 d half-life) has been used in sealed sources for 
industrial applications, for substitution of '^^Ir sources 
applied in conventional non-destructive testing (NDT) 
projectors, providing better quality of the radiographic 
pattern ||39l . The longer lived '^^ '^'^Lu are used as ra- 
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dioactive tracers in different processes BOl 1411 and the 
'^^Lu(''^Yb) decay is used for investigation ion implan- 
tation in semiconductors by using PAC technique i42ll . 
and for radiotracer studies of oil pipeline flow rates, re- 
finery column residence times, and the performance of 
a coal liquefaction pilot plant 1,43] . Deuteron induced 
reactions on ytterbium however, are not the most favor- 
able reactions for production of some of these radionu- 
clides as for instance ""^ Yb, '^^Lu and Yb can be pro- 
duced much more economically at reactors via (n,y ) re- 
actions. The optimal routes for the production of '^^Tm 
are the i*'^Er(p,n) |44| and i^^Er(d,2n) li23l . 
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Table 1: Decay characteristic of the investigated reaction products 



Nuclide 


Half-life 


E,(keV) 




Contributing reaction 


Q- value 
(keV) 


1''''SLu 


6.647(1 


112.9498 
208.3662 


6.17 
10.36 


"'^Yb(d,n) 
l^^Yb decay 


3959.3 


l'«Lu 


1.37 y 


78.63 

100.724 

272.105 


11.9 

5.24 
21.2 


'"Yli(d,n) 
^'^^Yb(d,2n) 
'"^"^YKd.Sn) 
'"^^^YKd.Sn) 


2689.9 
-3677.41 
-11142.05 
-23829.22 




6.70 d 


78.7426 
181.525 
810.064 
900.724 
912.079 
1093.63 


10.6 
20.6 
16.6 
29.8 
15.3 
63 


'^'Yb(d,n) 
"2Yb(d,2n) 
''''Yb(d,3n) 
''*Yli(d,4n) 
176YhfH fiti^ 


2493.6 

-5525.86 

-11893.17 

-19357.81 

-32044.98 


171(7. 


8.24 d 


667.422 
739.793 
780.711 
839.961 


11.1 
47.9 
4.37 
3.05 


^ '^Yb(d,n) 
''^'Yb(d.2n) 
'"^^YKd.Sn) 
'^-'Yb(d,4n) 
'"^"^YKd.Sn) 
^''^Yb(d,7n) 


2128.98 

^85.5 

-12504.96 

-18872.29 

-26336.93 

-39024.1 




2.012 d 


84.262 

193.13 

572.20 

985.10 

1054.28 

1138.65 

1280.25 

1341.20 

1364.60 


8.7 
2.07 
1.25 
5.4 
4.60 
3.49 
7.9 
3.15 
4.47 


^ "^Yb(d,2n) 
^^^Yb(d,3n) 
^^2Yb(d,4n) 
1^3Yb(d,5n) 
^^'*Yb(d,6n) 


-6465.7 
-13080.2 
-21099.6 
-27466.9 
-34931.6 




32.018 d 


109.77924 
130.52293 
177.21307 
197.95675 


17.39 
11.38 
22.28 
35.93 


'6^Yb(d.n) 
^'^'^Yb(d,3n) 
l'lYli(d,4n) 
l''^Yl)(d,5n) 
^^^Yb(d,6n) 


1567.07 

-13769.87 

-20384.37 

-28403.84 

-34771.16 




51.5 m 


178.87 
213.20 
239.22 
401.17 
1267.26 


2.5 
3.33 
7.7 
3.17 
3.87 


168Yb(d,3n) 
^™Yb(d,5n) 
'"Yli(d,6n) 
l'^2Yb(d,7n) 


-15151.7 
-30488.7 
-37103.2 
^5122.7 


l''''Yb 


1.911 h 


150.3 


20.5 


'^''YMd.p) 


3341.834 




4.185 d 


113.805 
282.522 
396.329 


3.87 
6.13 
13.2 


'^"^YKd.p) 
' ^6 Yb(d,p2ii) 
'"^^Tm decay 


3597.7842 
-9089.38 


169 


32.018 d 


109.77924 

130.52293 

177.21307 

197.95675 

307.52 

307.73586 


0.1739 
0.1138 
0.2228 
0.3593 
0.003 
0.1005 


'68Yb(d,p) 

170Yb(d.p2ii) 

''^'Yb(d,p3ii) 

''^2Yb(d,p4n) 

''^^Yb(d,p5n) 

^'^'^Yb(d,p6n) 

'*'lu decay 


4642.414 
-10694.52 
-17309.02 
-25328.49 
-31695.81 
-39160.45 


l'«Tm 


8.24 h 


398.9 
461.4 


87.9 
6.9 


'"Yl)(d,2p) 
'''''Yb(d,2pn) 
''^^Yb(d,2p3n) 


-2739.62 
-10204.25 
-22891.43 


l'>«Tm 


93.1 d 


79.804 
184.295 
198.251 
447.515 
720.392 
741.355 
815.989 
821.162 


10.8 
17.9 
53.8 
23.7 
12.0 
12.6 
50.3 
11.8 


""'Yb(d,2p) 

'™Yb(d,2p2n) 

''"Yl)(d,2p3n) 

l''2Yl)(d,2p4n) 

l''Yl)(d,2p5n) 


-1699.2 

-17036.15 

-23650.65 

-31670.12 

-38037.44 




9.25 d 


207.801 
531.54 


42 
1.61 


'68Yb(d,2pn) 
l™Yb(d,2p3n) 
l'^lYb(d,2p4n) 
l'^2Yb(d,2p5n) 
^^"^Yb decay 


-8539.87 
-23876.82 
-30491.32 
-38510.79 




30.06 h 


242.917 
297.369 
460.263 


35.5 
12.7 
4.12 


'68Yb(d,2p3n) 
^™Yb(d,2p5n) 
'^^Yli decay 


-24294.89 
-39631.84 



Abundance of isotopes in natural Yb (%): '''^Yb-O.lS, '™Yb -3.05, "'Yb -14.3, ™Yb -21.9, "^Yb -16.12, "''Yb-31.8, "*Yb -12.7 
The Q- values refer to formation of the ground state and are obtained from [18] 

When complex particles are emitted instead of individual protons and neutrons the Q-values have to be decreased by the respective binding 
energies of the compound particles: np-d, +2.2 MeV; 2np-t, +8.48 MeV; n2p-3He, +7.72 MeV; 2n2p-a, +28.30 MeV 
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Table 2: Measured cross-sections of the "o'Yb(d,xn)"'''i"'"2"'*.i7imj,i70,i69L„ reactions 



Energy ±AE (MeV) 


Cross-section +Acr (mb) 














39.1 


n, 1 




in.N 


1 














II : 




2 .■■ . S 


36.2 




1. ■■■.() 


]!.(' 




i.y 


.V. I./ 


12./ 






1 


1 1.2 


llN. \ 


21.1 


33.1 


0.8 


.SI. 8 


14.0 


246.4 


33.4 


482.5 


57.7 


459.6 


.52.7 


372.7 


41. H 


174.5 


20.3 


29.7 


0.9 


54.1 


15.5 


253.7 


31.2 


450.2 


51.2 


434.3 


49.8 


304.4 


34.2 


127.9 


15.0 


26.0 


1.2 


71.0 


24.9 


386.5 


45.3 


377.7 


43.0 


479.8 


54.7 


231.0 


25.9 


68.4 


8.4 


22.0 


1.3 


79.4 


23.0 


557.3 


65.1 


27S.3 


31.9 


395.5 


45.2 


196.8 


22.1 


29.7 


4.2 


17.2 


1.6 


182.7 


59.2 


377.1 


44.8 


253.5 


29.0 


259.4 


30.0 


82.4 


9.3 


6.7 


1.6 


11.4 


1.8 


286.3 


43.7 


98.9 


13.7 


106.1 


12.2 


80.4 


9.5 


11.8 


1.7 







Table 3: Measured cross-sections of the ""'Yb(d,x)''5 '^^Yb and "'"Yb(d,x)'"'"^'i**'i*^'^''5Tm reactions 



Energy ±AE (MeV) 


Cross-section +Acr (nib) 


1/3 Yb 


169Yb 








l63Tm 


39.2 


0.4 


72.5 


8.2 


371.4 


41.7 


1.3 


0.2 


12.5 


1.4 


14.3 


1.6 


2.1 


0.2 


36.2 


0.6 


78.8 


8.9 


305.4 


34.3 


1.3 


0.2 


12.4 


1.4 


8.7 


1.0 


1.1 


0.2 


33.1 


0.8 


75.8 


8.5 


236.0 


26.5 


1.2 


0.2 


16.0 


1.8 


6.6 


0.7 






29.7 


0.9 


70.3 


7.9 


190.7 


21.4 


1.5 


0.2 


12.5 


1.4 


6.9 


0.8 






26.0 


1.2 


65.5 


7.4 


100.5 


11.3 


0.9 


0.1 


8.6 


1.0 


5.3 


0.6 






22.0 


1.3 


61.4 


6.9 


44.0 


4.9 


0.9 


0.1 


6.8 


0.8 


4.3 


0.5 






17.2 


1.6 


65.9 


7.4 


10.0 


1.1 


0.6 


0.1 


5.8 


0.7 










11.4 


1.8 


74.2 


8.3 


1.5 


0.3 






1.9 


0.4 
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